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Tomonaga-Luttinger-liquid theory is formulated for metallic carbon nanotubes with open 
boundaries. Both cases of single- and multi-wall nanotubes are discussed. Based on this 
theory, spatial variation of the charge density from an edge is investigated with taking account 
of the shift of the chemical potential which expresses the carrier injection to the nanotube. 
The charge density has the spatially independent part and the oscillatory component. Roles of 
Coulomb interaction on the amplitude of the oscillation, the wavenumbers of it and the uniform 
component of the charge density are clarified. 
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§1. Introduction 

A carbon nanotube is composed of a coaxially rolled graphite sheetS and it's actual length is of the 
order of l^m or less than it. The material is characterized by two integers, (A^+, N^), corresponding 
to a wrapping vector along the waist, w = Nj^aj^ + N^a^, where a± = (iba/2, \/3a/2) are primitive 
lattice vectors of the graphite and |a-|-| = a is the lattice spacing. It has been shown that the 
carbon nanotubes have peculiar band structures.!' When — = mod 3, the metallic 
one-dimensional dispersions appear near the center of the bands. The low energy properties less 
than vq/R {vq : Fermi velocity, R : radius of the tube) are well described by taking into account 
only the metallic one-dimensional dispersions. Thus, the metallic carbon nanotubes are considered 
as the typical one-dimensional conductors. 

It has been well known that physical properties of the one-dimensional interacting electron sys- 
tems cannot be described by conventional Fermi-liquid-theory. Instead, the systems show the 
behaviors called as Tomonaga-Luttinger-liquid. The Tomonaga-Luttinger-liquid state is charac- 
terized by separation of the charge and spin degrees of freedom, and the anomalous exponents of 
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correlation functions, which depend on the interaction. Carbon nanotubes are one of the most 
promising candidates where such exotic correlation effects can be observed. The electronic states 
of the metallic carbon nanotube have been theoretically investigated with taking account of the 



perimentally, in the transport measurements for single- wall nanotubes (SWNTs), power-law 
dependences of the conductance as a function of temperature and of the differential conductance as a 
function of bias voltage have been observed in the metal-SWNT junctions and in the SWNT-SWNT 
junctions. These results have been interpreted to be due to tunneling between the Fermi-liquid and 
the Tomonaga-Luttinger-liquid in the former case, and that between the Tomonaga-Luttinger-liquid 
and the Tomonaga-Luttinger-liquid in the latter case. 

The correlation effects of the semi-infinitei'& and finite length carbon nanotubJl' have been 
studied theoretically. The local density of state of the semi-infinite SWNT calculated by the 
bosonization theory explains the transport experiments quantitatively. In ref.^, the local tunnel- 
ing density of states of the finite length carbon nanotube has been investigated at absolute zero 
temperature (T = 0). It has been found that the new energy scale, which reflects the spin-charge 
separation in the Tomonaga-Luttinger-liquid, appears in Coulomb blockade behavior in addition 
to usual charging energy and single-particle level spacing. However, in the above theories, the re- 
lationship between the original electron operator and the slowly varying Fermi field describing the 
low energy physics has not been clarified. The relationship is important for discussing the concrete 
physical quantities, especially the spatial dependence of those with rapid oscillation determined 
by Fermi wavelength. Note that the oscillatory component of the local density of states has been 
neglected in refs. ^ and |^ because the transport measurements are considered to observe the local 
density of states averaged over several lattice constants. In the present paper, we develop the 
bosonization theory of the metallic carbon nanotube with open boundaries based on the theory 
for one-dimensional systemB with paying attention to the relationship between the two kinds of 
electron operators. The theory is extended to the case of the multi-wall nanotube (MWNT). By 
utilizing the theory, the spatial variation of the charge density from an edge is investigated in the 
presence of the shift of the chemical potential expressing doping of carriers to the nanotube. It 
is found that the uniform component of the charge density and the wavenumber of the oscillation 
become smaller compared with the non-interacting case due to the effects of the interaction on 
the zero modes. The amplitude of the oscillation determined by the bosonic fluctuation is shown 
to be larger than that in the absence of the interaction. The amplitude for MWNT depends on 
the number of the metallic shells included in the MWNT, and the effect of the interaction on the 
amplitude vanishes when the number tends to infinity. 

The organization of this paper is as follows. In §.2, the bosonization theory with open boundaries 
is developed for the SWNT. The theory is extended to the case of the MWNT in §.3. The charge 



long-range 
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distribution of the SWNT and of the MWNT are calculated in §.4. Section 5 is devoted to summary 
and discussion. 

§2. Bosonization for SWNTs with open boundaries 

We formulate the bosonization theory of the metallic SWNT with open boundaries. As a model 
of the metallic carbon nanotube, we consider the {N, N) armchair nanotube. 

2.1 Non-interacting case 

We consider the armchair carbon nanotube with the length, L, and the radius, R = \/?>Na/ {2it), 
schematically shown in Fig. 1. The Hamiltonian in the tight-binding model is written as 7i = 
+ "^int, where TL]^ is the kinetic part, 

= -t 51 + - + a-p^s{r - pa_)|, (1) 

r,s,p 

and 7i\nt describes the mutual interaction, 

Wint = \Y.Y.Y.U{r-r' - {p-p')d/2)al,{r)al,^Ar')a,,,Ar')apAr)- (2) 
r,r' s,s' p,p' 

Here t denotes the hopping integral between the nearest-neighbor atoms, and ajj^^{r) is the creation 
operator of the electron with spin s = it at the location r — pd/2 where r = (alj^/Sam) or 
{a{l + l/2),^/3a{m -\- 1/2)) (/, m : integer), p = it and d = (0,a/-v/3). The interaction, U{r), is 
given by U{r) = e^/ I^k^ Cq + x"^ + iR"^ sin^(y/2i?)|i) with r = {x,y) where e, n and oq are the 
electric charge, dielectric constant and short-range cut-off of the interaction, respectively. 

At first, we consider the non-interacting case. In order to diagonalize T^k, we solve Schrodinger 
equation, 

-t{f^'{r) + f^^{r-pa+) + f^^{r-pa^)} = E^{K)f^{r), (3) 

with the boundary condition, 0^(0, y) = (l/^{L,y) = and (/'^(^•) = (p^-^if + '^^) where w = 
A^(a+ -|- a_) = (0, \/3Na). The solution and eigenvalue are obtained as follows, 

'^k(^) = ]f^A-'P{K)smK,xe'^yy, (4) 
E^{K) = -rt\aiK)\, (5) 

where r = ±, A'"'+{K) = 1/^/2 and ^^"(K) = ra{K) / {V2\a{K)\) with a{K) = 1 + 
2 cos(-fCa;a/2)e'^2/v^"/2_ The wavenumbers, and Ky, are quantized as Kx = UxTt/L (ux = 
1,2, - ■■, L/a) and Ky = 2imy / {\/?>N a) {uy = 1, 2, • • • , 2A^). We note that the energy disper- 
sions with Ky = 271 /{VSa) = K^ touch the Fermi energy {E''{K) = 0) at Kx = 27r/(3a) = 
Kq, and the others do not cross the Fermi energy. The Hamiltonian is written as Wk = 
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Fig. 1. Carbon atoms in the armchair nanotubes with length, L, and the radius, R, where the x(y) axis is along 
the tube (waist). Here a± are two primitive lattice vectors of graphite, \a±\ = a, and d — (0, a/\/3). The hexagon 
delineated by the thick line is the unit cell and the black (white) circle denotes the sublattice p — +(— ). 



J2k ^r,s EriK)cl^^{K)cr,s{K) where ap^s{r) = Y^K^r 't)''j^{'f')(^r,s{^)- When the chemical po- 
tential, /i, stays the energy bands with Ky = Ky (see Fig. 2), the charge density in the absence of 
the interaction is given at T = for the semi-infinite case, L ^ oo, as 

■r-^ t / N (K^ — K^)a a (sm2K:X — sm2K_x) 
{: T.<s{r)a,Ar) :) = , (6) 

where K± ~ Kq it fi/vo and vq = \/^ta/2. The charge distribution has two parts. One is the 
uniform part due to the shift of the chemical potential and the other expresses the oscillation 
(Priedel oscillation). We find that there exit the two kinds of oscillation with the wavenumbers, 
2K±, and those are out of phase. The oscillatory behavior is absent for fi = 0. 

We apply the bosonization procedure on the system with the open boundarieJil' to the present 
model. Energy dispersion is linearized by substituting K = (Kq + k,Ky) into E^'{K). When we 
define the right/left moving electron operator, a^/2,,s(-^o + k), as 

f c_ siK), k> 
\ c+,,(K), A:<0 

aL,s{Ko + k) = i (8) 
[ c^,siK), k<0, 

the Hamiltonian, "H^, and the electron operator are respectively written as follows, 

= 51^0^ {(^R,s(Ko + k)aR4Ko + k)- a^^Ko + k)aL,s{Ko + k)] , (9) 
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K_ 1 
K^(n/a) 



Fig. 2. Energy dispersion, E{K), of the (5,5) armchair nanotube as a function of Kx- Fermi energy in the absence 
of the chemical potential, n, is E(K) = and K± ~ ifo ± l^-fvo where Ko = 27r/(3a) and vo = V3ta/2. 



ap,s{r) = V ^ E sin(^o + A;)xe'^^ {paR,s{Ko + k) + ulA^o + k)} . (10) 

k 

Here we define the slowly varying Fermi field, '(pR/L,a,s{^)j ^ 

(11) 
(12) 

(13) 

(14) 

where a = it is the index corresponding to valleys for the case for periodic boundary condition. 
The electron operator is expressed by the slowly varying fields as 





A^) = - 


—1 

\/2L 


Y^e^'-aRAKo + k), 

k 




A^) = ■ 


—i 
\/2L 


Y,^^''aLAKo + k), 

k 


i'L- 


A^) = - 


i 

\/2L 


Y^e-'^^^'aRAKo + k), 

k 


V'R,- 


A^) = - 


i 

v/2L 


Y^e-'^^'^aLAKo + k), 



«^''^(^) = V {^'"""^ Hi^,+,.(x) + V'L,+,.(x)] + e-'^°- [pVl,-,.(x) + ^R,-A^)]} . (15) 

The expression is the same as that in case of the periodic boundary condition. The slowly varying 
field operators are not independent, but satisfy, 

tpL,+A^) = -4^R-,s{-x), (16) 

^L-A^) = -i^R,+A-x)- (17) 
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Therefore we can actually work with the right moving operator only. The boundary condition, 
Q^p,s(0, y) = 0, is automatically satisfied due to eqs. (16) and (17). However, the condition, 
ap^s{L,y) = 0, implies, 

V'i?,a,s(-i) =e'2"^°^^R,„,.(L). (18) 
So we can regard the field ipR,a,s{^) as defined for < x < 2L with the boundary condition, 

^R,o,,s{x + 2L) = e-*2"^«^Vi?,a,.(x) 

^gi2aW3^^_^^^(^)^ (19) 

where = 0, ±1 is defined as L = (Sn/, + v)a (hl '■ integer satisfying L > 0). In terms of the right 
moving operators, eqs. (9) and (10) are expressed as 

'H], = voJ2 [ dxV'Jj_^^^(-i5a;)V'i?,«,s, (20) 
ap,s{r) = e'^^X^^p(i+°)/2 {ei-^»^V'fl,„,,(x) - e--^°-Vfl,a,.(-^)} ■ (21) 

The above Hamiltonian and the right moving electron operators are straightforwardly bosonized 
by utilizing Haldane prescription.© The bosonized form for 'ipR^a,s is given by 

v27ra 

where AA^^^^ is the extra electron with the index (a,s) and satisfies, [Oa^g^ ^N^i ^^i] = \6a^a'^s,s'- 
The operator, rja^g, is Majorana Fermion satisfying {?/o,s, r/Q.',s'} = '^.da^a'^s^s' and a is the ultra 
violet cut-off of the order of 1/R. The function, <f)ct,s{x)-, is given by 

^o.Ax) = E , (23) 

and satisfy 4>a,s{x + 2L) = 4>a,s{x) because q = rm/L (n : integer). Here ba,s{Q) is the bo- 
son operator, so satisfies [^o,s('7); ^j^/ g'('7')] ~ ^a,a'^s,s'^q,q'- The density operator, PR^a 
V'lj,a,s(^)V'i?,a,s(2;), is given by 

A7V„,^ + az^/3 1 

PR,aAx) = + ^dx<t)a,s{x), (24) 

and pL,a,s{x) = pR-a,s{—x). The bosonized form of is given as 

a,s 

= ^ E(^^".« + f + E E ^^096^(9) (25) 

Before including the mutual interaction, we change the variables of (a, s) to those for charge/spin 
(i = p/^) degree of freedom of symmetric/antisymmetric {6 = +/—) combination between the 
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valley index. We define 



a,s 


(26) 


bp-{q) = ^Y^^^^Aq): 


(27) 


ba+{q) = ^^s6„,s(g), 


(28) 


\ 

ba-{q) = i^'Yasba,s{q), 


(29) 


ANp+ = J2 ^Na,s, 

a,s 


(30) 


a,s 


(31) 


a,s 


(32) 


a,s 


(33) 


, [9a,s,ANp+] = i, [0«,„AiVp_] 


= ia, [6'a,s, A7Vcr+] = is and 



[6a,s, AA^'o--] = ias. Since the eigenvalue of ANa^s is an integer, that of ANjs, Qjs has a constraint, 



Qp+ + aQp- + sQcr+ + asQcr- = 4 x integer. 
In terms of the new variables, and i^R,a,s{x) are respectively expressed as follows, 



(34) 



j,5 j,5 q>0 

= E 

r]a,s ^-ie^.s^i^^ANo++a(ANo-+iu/3)+sAN^++asAN^-} 
X ^h^'i'P^ {x)+a(j>p- (a:)+s0a+ {x)+as<j>a- (x)} 



with 



PjS 



» = T.\[^{^'"'%s{q)+h.C.}. 



(35) 



(36) 



(37) 



2. 2 Effects of interaction 

We take into account of the mutual interaction, eq.(2). Among the various interaction processes, 
the term with the strongest amplitude is written in terms of slowly varying Fermi fields as J 

V{0) 



[ dx \^[pR^a,s{x) + PL,a,six)]\ 
^0 [^s J 



ViO) 



dx 



+ E E PR,ci,s{x)PR,a',s'{- 



(38) 



a,s a's' 
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where V{0) = {2e'^ / k) ln{Rg / R) where Rg ^ L characterizes the large distance cut-off of the 
Coulomb interaction. The other interaction processes, whose couplings scale as a/(27ri?) and are 
extremely small, have been known to play crucial roles in the absence of the carrier doping, i.e., 
at half-filling.i'i'& Since we discuss the case for away from half-filling by introducing the shift of 
the chemical potential, those are safely neglected. The above interaction has a following bosonized 
form, 



q>0 q>0 



(39) 



Since the interaction term, eq.(39), is expressed by {p+) mode only, the Hamiltonian except {p+) 
have diagonalized form, Ti-^s- "^^^ Hamiltonian of {p+) mode, Ti.p+, is diagonalized by Bogoliubov 

a result, Ti.p+ are given as 



transformation, bp^{q) cosh (pbp^{q) — sinh ipb^^{q) with e^'^ = 1/^1 + AV{0)/{ttvo) = Kpj^. As 



Wp+ = ^^(AiV,+ )2 + Y vp+qbl^iq)bp+{q), (40) 

where Vp^ = v^/Kp^ and fp+A? = Vq/K^^. The quantity, (j)p^{x), is transformed as 

(j)p+ (x) cosh ip(j)p+ {x) - sinh ip(j)p+{-x) , (41) 

in eqs.(36). Prom eq.(40), the charge susceptibility per unit length is easily derived as 4/{nVp^N)- 
The quantity is suppressed by the long-range Coulomb interaction because Kp^ < 1. For typical 
metallic nanotubes, the value of i^'p+ is estimated as about O.2JI' which leads to strong suppression 
of the charge susceptibility as A/^ttvq) x 0.04. 

§3. Extension to multi-wall nanotubes 

In this section, we extend the bosonization theory with open boundaries on the SWNT developed 
in §.2 to the case of MWNT.© We consider the MWNT where Nm metallic graphite shells with 
radii Ri < R2 < ■ ■ ■ < Rnm included. For simplicity, we consider the case where all the metallic 
shells consist of armchair nanotubes. The insulating shells in the MWNT can be incorporated in 
space-dependent dielectric constant. The electron operator for the n-th shell, ap^s^„(r), is given as 
follows. 



where 27ri?n = \/^Nna and 



^p(i+")/2 {e-^o-^H,„,,,„(x) - e-^-^»^^R,„,,,„(-x)} , (42) 



(x) = p-iec,s,npi^{ANp+„+a{ANp-„+iiy/3)+sANa+„+asAN„-„} 

^g^{</ip+„(x)+a</ip_„(a;)+S(/)CT+n(a;)+os<iiCT_„(a;)} ^-^g^ 
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where {'na,s,n,Va' ,s'n'} = '^Sa,a'^s,s'Sn,n', « ~ ^/^Nm, [da,s,n, '^^p+n'] = ^Sn,n' , [(^a,s,n, '^Np^n'] = 
iOf^n.n', [(^a,s,n, ^N„+n'] = is5n,n' , [9a,s,n, ^N^_n'] = ^as5n,n' , and 

^M^) = E a/^ U''''-'^''%5n{q) + h.c.) , (44) 
with [bjSn{Q),b^j's'n'(^')^ ~ ^j,j'^s,S'^n,n'Sq,q' ■ The Hamiltoiiian is written in terms of the bosonic 



variables as 



= 1^ E E(^^.5n + ihA,-f + E E E ^og6j,„(g)6,.„(g), (45) 

n=l j,5 n=l j,5 q>0 

"^int = E '^X^^ {pR^a,s,ni^)PR,a',s',m{x) + PR^a,s,n{^)PR,a',s',m{-x)} 

n,m=l ■'^^ a,s a',s' 

TTVo 

= -^J^ E ^nm^N p+n^ p+m 
n,m=l 



E E 1 bl+n{(l)bp+m{q) - 2 (W(9)&p+m(9) + b\+m{(l)bl+n{l)) \ ■ (46) 



g>0 n,m=l 

Here C/nm = 4V„m/(7'"Wo) = {8e^/(7rt>oKnm)} ln(-Rs/^nm) expresses the interaction between the n- 
th and m-th shell with Knm and Rnm being the dielectric constant between the n-th and m-th shell 
and the "mean radius" of i?„ and Rm introduced in ref.^, respectively. The modes except {p+) 
are already diagonalized. By using the orthogonal matrix satisfying X]m=i Unm^m,j = ^njdj, the 
Hamiltonian for (/>+) mode, 'Hp+ is written as 

^p+ = |^E(l+*)(A^P+.)' 

Nm f > 

+ E^ogE (1 + ^iH+Mbp+M - ^i{bp+,iQ)bp+M + bUMbUM) ' (^7) 
g>o j=i L ^ 4 J 

where AiVp+„ = J2f=i'^nj^Np+j, bp+n{q) = T,f=i'^njbp+j{q) and gj is the eigenvalue of Unm- 
Equation (47) is diagonalized by Bogoliubov transformation, bp+j{q) —>■ cosh (pjbj{q) — sinhipjbj{q) 



as 



= E I ^(A^P+.)' + E mb]{Q)bj{q) \ , (48) 

j = l [ q>0 



where vj = vq/Kj and WjAr = vq/KJ with i^j- = 1/ ^yi + gj. In terms of 6j(g) and bj{q), (j)p+n is 
expressed as 

Nm 

4>p+nix) = ^ Tnj {coshv3j(/)j(x) - sinh (^j</)j (-x)} , (49) 
'^.•(^) = E + /^-c-} • (50) 
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§4. Charge distribution 

Based on the bosonization theory formulated above, the distribution of the charge density from 
the one boundary is discussed with taking account of the shift of the chemical potential, /x. The 
roles of the zero modes and of the bosonic long wavelength fluctuations are clarified. 

At first, we consider the SWNT. The term expressing the shift of the chemical potential, — /xAiVp+, 
is added to the Hamiltonian, eq.(40). Using eq.(21), the charge density is given as follows, 



AN 



The first and the second terms are spatially independent and given by 



(52) 



s,a 



On the other hand, the third and fourth terms are given by 



Y: {e-^^"^"-(<„,,(x)V«,a,.(-x)) + e^^"^o-(4,„,,(-x)V^^,„,,(x))} 

a,s 

1 {Kp++3)/8 



27ra 

oo 

xn 



sinh^(7ra/2L) 



sinh"'(7ra/2L) + sin^(7rx/L) 



sin^(7rx/L) 
sinh2(/7r2^p+/L) 



1 + 



sin^(7rx/L) 
sinh2(/7r2^o/i^) 



n-3/4 



(53) 



where f{2x) = tan"^ {sin(27rx/L)/(e''^/^ - cos 27ra;/L)}, = Vp+/{2TrT) and = vo/{2TrT). 
Here (• • •) in the r.h.s. in eq. (53) expresses the average in terms of the zero modes. Since the first 
line in the r.h.s in eq.(53) expresses the oscillation of the charge density, the zero modes determine 
the uniform shift of the charge density by the chemical potential, eq.(52), and the wavenumbers of 
the charge density oscillation. On the other hand, the second and third lines are due to the bosonic 
fluctuation and express the amplitude of the oscillation. 

We consider the semi-infinite case, L ^ oo. In this case, the average in terms of zero-modes is 
easily calculated, 

4L/i 



(AiVp+) 



TTVp+N 



^^±i^{ANp++a{ANp_+4.u/3)+sAN^++asANa-}\^ ^ ^±i2nx/vp+N 

and f{2x) = (7r/2)sgn(x). As a result, eq.(51) is calculated as follows 



(54) 
(55) 



a 1 



TvNvp^ N N 27ra 



sin(2Ko + 



2/x 

Vp+N' 



)x — sin(2Ko 



2^ 

Vp+N' 



■)x\Asix), 
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(56) 

Note that eq.(56) together with eq.(57) in case of T = and Kp+ = 1 reduces to the non-interacting 
case, eq.(6). The deviation of the uniform charge density, fia/{TTNvp^]\f) = ij.aKp_^ / {nvo) , and the 
shift of the wavenumber of the oscillation, Iji/vpj^N = '^jJ.Kpj^/vo, are smaller than those in the 
absence of the interaction. Both quantities are about 0.04 as large as those in the absence of the 
interaction for the typical SWNT. The fact means that Coulomb interaction prevents the carriers 
being injected into the nanotube, and has the same origin as suppression of the uniform charge 
susceptibility. The amplitude of the oscillation, As{x), is large and decays slowly compared to 
the non-interacting case as is shown in Fig.3. Here, the amphtude of the oscillation is shown 
as a function of x/a for T/Ec = (solid curves) and T/Ec = 0.04 (dotted curves) with Ec = 
vo/a. For each temperature, the upper and lower curves express the amplitude in case of Kp^ = 
0.2 and that in the absence of the interaction, i.e., Kp^ = 1, respectively. The temperature, 
T/Ec = 0.04, corresponds to room temperature for N = 10. Note that the amplitude delays as 
exp {-{K^+ + 3)x/(4Co)} for x » Cp+- 




Fig. 3. The amplitude of Priedel oscillation of SWNT, As{x), defined in eq.(57) as a function oi x/a for T/Ec = 
(solid curves) and T/Ec = 0.04 (dotted curves) where Ec = vo/a. For each temperature, the upper (lower) curve 
expresses the amplitude in case of Kp+ = 0.2 {Kp+ = 1). 

Next we consider the case of the MWNT. The charge density of the n-th shell is obtained for 
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L — oo as, 

(: I^«p,s,n('^)op,5,n(r') :) 



s 



-^NnJ-i^ VjN 



X 



where = Vj/{2TrT). Because of the weak logarithmic dependence of Unm on Rnm, we can obtain 
sensible results by the approximation Unm = U. The approximation leads to gi = NmU, gj = 



for j = 2 • • • Nm and Fni = l/V-^M- As a result, the charge density is obtained as 

{■■ E4,«,n(^)«P,.,n(r) :) = - -^-^ |sin(2i^o + - sin(2i^o - ^)x] Am(x), 

(59) 

, , , / ~a\ iKi-i)/mM)+i r sinh(x/ei) 1 -^1/(4^^) r sinh(x/eo) 1 -i+V(4A^m) 

where i^i = l/Vl + NmU, = v\/{2'kT) with = vq/Ki and t^ijv = vq/K^. The amplitude of 
the oscillation, ^^(ic), becomes small with increasing Nm as is shown in Fig.3, where the quantity, 
Am{x), is shown as a function of x/a for T/Ec = (solid curves) and r/£?c = 0-04 (dotted curves) 
with Rnm = 6nm and Knm = 1-4 for several choices of Nm- For each temperature, the curve 
from top to bottom corresponds to the case of Nm = 1,2, 5, oo. We note that the effects of the 
interaction on Am{x) disappear in the limit, Nm — ^ oo. 

§5. Summary and discussion 

In the present paper, the bosonization theory with open boundaries was developed for both 
SWNT and MWNT. We payed attention to the relationship between the electron operator of the 
tight-binding model and the slowly varying field. Based on this bosonization theory, we investigated 
the distribution of the charge density from an edge with taking account of the shift of the chemical 
potential expressing carrier doping to the nanotube. 

The bosonized Hamiltonian is written by the sum of symmetric/antisymmetric combination of 
the charge/spin excitation. Each mode has two contribution, zero modes and bosonic fluctuation. 
Effects of the Coulomb interaction with the strongest amplitude appear in only the symmetric 
charge excitation. These facts are the same as the case of the periodic boundary condition. However, 
the Hamiltonian of zero modes are expressed in terms of only the extra number. This fact is due to 
the open boundary condition, and different from the periodic boundary condition, where the zero 
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Fig. 4. The amplitude of Friedel oscillation of MWNT as a function of xja for T/Ec = (solid curves) and 
T/Ec = 0.04 (dotted curves) with Rnm ~ 6nm and Unm = 1.4 for several choices of Na-i. For each temperature, 
the curve from top to bottom corresponds to the case of Nm — 1, 2, 5, oo. 



mode Hamiltonian is written by the extra number and the current. The charge distribution in the 
presence of the shift of the chemical potential has two component. One is spatially independent and 
the other shows oscillatory behavior which vanishes in the absence of the chemical potential. The 
zero modes determine the magnitude of the spatially independent component and the wavenumbers 
of the oscillation. The long-range Coulomb interaction strongly suppress the both quantities. This 
means the fact that the interaction prevents carriers from being injected into the nanotubes. For the 
typical SWNT with Kpj^^ = 0.2, the both quantities are about 0.04 as large as those in the absence 
of the interaction. The quantity, the shift of the chemical potential, corresponds to the applied gate 
voltage or difference between the work function of the nanotube and that of the substrate material. 
The carrier doping from the gold substrate has been observed by scanning tunneling spectroscopy 
(STS)0'0) and scanning tunneling microscopy (STM).lll' In the STS experiments, the asymmetry 
of the density of states has been observed, and it is maintained to be due to the difference between 
the work function of the nanotube and that of the gold(lll) substrate. In the present formalism, 
the shift of the density of states is just given by the chemical potential even in the presence of 
the Coulomb interaction, which is the same as the conclusion in refs.O and Q8l However, the 
deviation of the Fermi wavenumber, 6k, in the presence of the chemical potential is given by not 
the simple form of ^/vq but fi/vp^N' = fJ-^^p+f'^o as is seen in eq.(56). On the other hand, the STM 
experiment observes 5k ~ h/vq. The discrepancy seems to need the further theoretical study on the 
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Coulomb interaction of the nanotube. The amphtude of the oscillation determined by the bosonic 
long wavelength fluctuation shows the non-integer power law behavior as a function of the spatial 
coordinate at T = 0. We found that the amplitude is enhanced by the interaction and it becomes 
smaller with increasing the number of the metallic shells, A^m, in the MWNT and identical with 
that in the absence of the interaction in the case, Nm oo. Finally, we note that the present 
result of the charge distribution are also valid for the strong impurity potential. 
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